March 1993 



Sterile Neutrinos as Dark Matter 



Scott Dodelson 1 ' and Lawrence M. Widrow 2 ' 

1 NASA/Fermilab Astrophysics Center 
Fermi National Accelerator Laboratory 
P.O. Box 500, Batavia, IL 60510 

2 Department of Physics, Queen's University 
Kingston, Canada, K7L 3N6 
and 

Canadian Institute for Theoretical Astrophysics 
University of Toronto, Toronto, Canada, M5T 1A7 

ABSTRACT 

The simplest model that can accomodate a viable nonbaryonic dark matter 
candidate is the standard electroweak theory with the addition of right-handed or 
sterile neutrinos. This model has been studied extensively in the context of the 
hot dark matter scenario. We reexamine this model and find that hot, warm, and 
cold dark matter are all possibilities. We focus on the case where sterile neutrinos 
are the dark matter. Since their only direct coupling is to left-handed or active 
neutrinos, the most efficient production mechanism is via neutrino oscillations. If 
the production rate is always less than the expansion rate, then these neutrinos 
will never be in thermal equilibrium. However, they may still play a significant role 
in the dynamics of the Universe and possibly provide the missing mass necessary 
for closure. We consider a single generation of neutrino fields (u^, vr) with a 
Dirac mass, /i, and a Majorana mass for the right-handed components only, M. 
For M > /i we show that the number density of sterile neutrinos is proportional 
to /j?/M so that the energy density today is independent of M. However M is 
crucial in determining the large scale structure of the Universe. In particular, 
M ~ 0.1 — 1.0 keV leads to warm dark matter and a structure formation scenario 
that may have some advantages over both the standard hot and cold dark matter 
scenarios. 
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The COBE DMR experiment's recent detection of large-scale anisotropy in 
the cosmic microwave background 1 has considerably strengthened the view that 
the large scale structures seen today evolved from very small primeval density 
inhomogeneities. Still, the two primary ingredients which dictate how structure 
forms, namely the nature of dark matter and the shape of the primeval fluctuation 
spectrum, remain unknown. 

The best studied and perhaps most successful model for structure formation 
is known as the cold dark matter (CDM) theory 2 . In the standard CDM model, 
the Universe is assumed to be spatially flat (Q = 1) with 90 — 95% of the mass 
density in dark matter and the balance in baryons (5 — 10%) and photons and 
light neutrinos (^C 1%). Primeval fluctuations are generated during inflation and 
are Gaussian with a scale-invariant spectrum. CDM, with the additional assump- 
tion that galaxy formation is 'biased' to occur first at the highest peaks in the 
density fluctuation spectrum can successfully explain galaxy-galaxy and cluster- 
cluster correlation functions on scales of order 1 — 5 Mpc and is at least consistent 
with the morphology of galaxies. However, CDM now appears to be inconsistent 
with various sets of observational data. Perhaps its greatest difficulties come with 
large scale structure data such as the APM galaxy survey 3 , which suggest more 
power on large scales than standard CDM model predictions. On small scales, 
the observed pairwise velocity dispersion for galaxies appears to be significantly 
smaller than those predicted by CDM 4 . 

One alternative 5 which has recently received a fair bit of attention is cold + 
hot dark matter (C+HDM). HDM is usually taken to be a light neutrino with 
m v = (92fij,/i 2 ) eV where H = lOO/i km/sec/Mpc is the Hubble parameter. In 
models with HDM alone, the processed fluctuation spectrum is characterized by 
the typical distance a neutrino travels over the history of the Universe, \ v ~ 
40 (30 eV /m v ) Mpc. This in turn sets the mass scale below which damping occurs 
due to free-streaming, M FS = Anp (A^/2) 3 /3 ~ 3 x 10 15 (30 eV/m^) 2 Q- 1 M . 
In HDM models, the first structures to form are pancake-shaped objects of size 
\ v with smaller scale structures such as galaxies and clusters forming later via 
fragmentation. However, we know from the galaxy correlation function, that the 
scale which is just becoming nonlinear today is around 5/i _1 Mpc. Essentially, the 
problem with HDM alone is that is too large: If galaxy formation occurs early 
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enough to be consistent with high-redshift galaxies and quasars, then structures on 
5/i _1 Mpc will overdevelop. The hope is that C+HDM will combine the successes of 
both models. In fact, a survey 6 of models with various amounts of hot dark matter, 
cold dark matter and baryons points to ^baryon — 0.1, ^CDM — 0.6, Vt v = 0.3 and 
a Hubble constant of h = 0.5 as the best fit model for microwave anisotropy data, 
large scale structure surveys, and measures of the bulk flow within a few hundred 
megaparsecs. 

As appealing as C+HDM may be for large scale structure phenomenology, it 
is somewhat unpalatable from the point of view of particle physics. Since there 
are no stable, neutral, massive particles in the 'standard model' for electroweak 
interactions, the existence of nonbaryonic dark matter implies new physics. Given 
that the existence of the baryon-antibaryon asymmetry also requires new (and 
probably distinct) physics, it seems already a great coincidence that O^M and 
^baryon be as close as they are 7 . Two types of dark matter imply further additions 
to the standard model with yet another coincidence in order to have ^hdm> ^cdm> 
and ^baryon an within one or two orders of magnitude of each other 8 . 

By far the simplest dark matter candidate, at least from the point of view of 
particle physics is the neutrino. Massive neutrinos require only the addition of 
right-handed or sterile neutrino fields to the standard model. In fact, it is the 
absence of right-handed neutrinos that seems contrived in light of the fact that all 
other fermions in the standard model have both left and right-handed components. 

Here we focus on the possibility that sterile neutrinos are the dark matter and 
that they are somewhat heavier but less abundant than the usual HDM neutrino. 
Such a 'warm' dark matter particle may have advantages for structure formation 
over both hot and cold dark matter scenarios. Our work is similar in some respects 
to that of Bond, Szalay, and Turner 9 who consider a particle that is in thermo- 
dynamic equilibrium at early times but decouples before ordinary neutrinos do so 
that <?*, the number of effectively massless degrees of freedom, is relatively high 
(g* > 100). Warm dark matter has been for the most part been ignored, to a 
large extent because there have been no compelling candidates proposed thusfar. 
In part, the motivation for this work is to propose a 'realistic' warm dark matter 
candidate. 

For simplicity, we consider only one generation of neutrinos. The mass terms 



3 



for the neutrinos are then 10 : 

C = n (J^J v l vr + Mu R u R + h.c. (1) 

where is the standard model Higgs field with ((f)) = v. The usual HDM case, 
wherein the active neutrinos constitute the dark matter, corresponds to {// = 92h 2 eV, M <C 
or {/i 2 /M = 92h 2 eV, M 3> /x}. When sterile neutrinos are the dark matter, the 
relevant mass is M. At tree-level, vr couples only to vl and therefore the most 
efficient way to produce sterile neutrinos 11 ' 12 ' 13 is via oscillations vl — > vr. The 
probability of observing a right-handed neutrino after a time t given that one starts 
with a pure monoenergetic left-handed neutrino is sin 2 26 m sin 2 vt/L where 6 m is 
the 'mixing angle', L is the oscillation length, and v is the velocity of the neutrinos. 
In vacuum, and with fx -C M (see-saw model) 6m = f^/M and L = AE/ (M 2 — /i 2 ) 
where E is the energy of the neutrinos. In the early Universe, the observation 
time t is replaced by the interaction time for the left-handed neutrinos. Recent 
work 14 ' 15 ' 16 has fine-tuned this picture taking into account the effect of finite den- 
sity and temperature on the mixing angle. 

Here we are interested in the case where the right-handed neutrinos are pro- 
duced at temperatures of order 100 MeV though the production rate is never so 
fast that they equilibrate. We begin with the Boltzmann equation for the sterile 
neutrinos: 



d d 
J± _ HE^- } f s (E,t) 
dt dE ' JbK " ' 



^sin 2 (2^ M (E,t)) F(E,t) 



f A (E,t) (2) 



where fs and Ja are the distribution functions of the sterile and active neutrinos. 
In the epoch under consideration (T>1 MeV) the left-handed neutrinos are in 

thermal equilibrium so that Ja = [e E ^ T + 1 j — (e p / T + 1 j . The quantity 
in square brackets is the probability per time of an active neutrino converting 
into a sterile one 16 where we have used the fact that for parameters of interest, 
the collision time is always much greater than the oscillation time (i.e. sin 2 vt/L 
averages to 1/2). The mixing angle and the collision rate are 17 

W(2^) ^ 2 + |(crg/ ^ ) + (M/2)|2 ; (3) 

where c ~ 4 sin 2 (2%) /15a ~ 26. 
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To get a feel for when and how many sterile neutrinos are produced, we derive 
the equation for r = ngjn^ where rii = 2 J cfipfi/ (27r) 3 is the number density of 
sterile (active) neutrinos with % — s (i — A). Changing the time variable from t to 
a, the Robert son- Walker scale factor and integrating Eq. (2) over momenta, one 
finds that 

dr 7 d\ng* 

= h r (4) 

din a if din a 

where 

7 "^/|^ sin22%(p ' r)r(p ' r) 57^TT' (5) 

and we have used the fact that g*a 3 T 3 = constant. For constant, j/H gives the 
number of sterile neutrinos, relative to the number of active neutrinos, that are pro- 
duced in each log-interval of T. Substituting Eq. (3), using H = 1. 66 gl' ' 2 T 2 /VftpianckJ 
and taking the limit M ^> /i, we find that 



jL 
H 



13 

~T/2 

g* 



y dy 

(ev + 1) (1 + x 2 y 2 Y 



(6) 



where x = 78(T/GeV) (keV/M). Taking g* = 10.8 and doing the integral nu- 
merically, we find that j/H reaches a peak value of 1.9 (/x/eV) (keV/M) when 
x ~ 0.19 or T = T max ~ 133 (M/keV) 1/3 MeV and falls off as T 3 for T < T max 
and T -9 for T ^> T max . Evidently, the number density in sterile neutrinos is 
proportional to M _1 so that the energy density is independent of M. Note also 
that most of the neutrinos are produced when the Universe has a temperature 
T ~ T max . As will be discussed below, our calculations simplify if we can assume 
that g* is constant. Since g* changes abruptly at T ~ 200 MeV and varies slowly 
for 200 MeV > T > 20 MeV, this assumption will be pretty good for M < keV 
but breakdown for masses much larger than this. 

Our interest is in the structures which form in a z/R-dominated Universe and 
we therefore require the full sterile neutrino distribution function. Here, we make 
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the assumption that g % is constant. Using dfs/dt = 



-HTdfs/dT and the identity 



T 




E 



+ E 




T 



T 



( 



dfs 



or 



) 



E/T 



(7) 



and changing the integration variable from T to x one finds 



oo 



fs 
Ja 



7.7 



1/2 




(1 + yV 2 ) 2 



dx' 



(8) 



X 



where y = S/T. In general, the right hand side of Eq. (8) is a complicated function 
of E and therefore will have a different energy dependence than ] a- There is no 
reason to expect otherwise: high energy and low energy neutrinos oscillate at 
different rates. Moreover, these rates change with temperature. HOWEVER, for 
T <C T max the lower limit of the integral can be set to zero and the right hand side 
of (8) becomes independent of E and T. In this limit, the integral is easily done 
and we find 



fs has the same functional form as Ja and therefore VLs/^v — [Mjm v ) (fs/fA)- 
From the relation m v j^l v ~ 92/i 2 eV we find that Vis = 1 for /i = 0.22/i eV 
where we have again set g* = 10.8. Finally, we note that the contribution of 
sterile neutrinos to the energy density of the Universe at the time of primordial 
nucleosynthesis 18 must be < 0.5 times the contribution of a light neutrino species 
if standard big bang nucleosynthesis 19 is to be believed. This in turn implies that 
M > 200/i 2 eV; that is, if sterile neutrinos are the dark matter then they are 
necessarily more massive than the standard HDM. 

How do perturbations evolve when a sterile neutrino species is the dark matter? 
Several guiding principles help us understand the processed power spectrum. First, 
structure within the horizon grows only after the dominant component of matter 
becomes nonrelativistic and therefore the size of the horizon at matter-radiation 
equality A#(a = a eq ) = a eq Jq EI1 dt'/a(t'), defines a characteristic scale. Second, 
perturbations on scales smaller than the Jeans length Aj = (^^pianck// 5 ) 1 ^ 2 
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(where v s is the speed of sound) oscillate like pressure waves. Finally, for neu- 
trinos, or any particle which is not completely non-relativistic, perturbations on 



exponentially damped. With the distribution function in Eq. (9), one can cal- 
culate these scales for sterile neutrinos. Figure 1 shows the relevant mass scales 
(= 47rp(A/2) 3 /3) as a function of the scale factor for the sterile neutrinos dis- 
cussed here and for an ordinary light neutrino dark matter candidate. For light 
neutrinos, the damping scale and the horizon scale at equality are roughly equal 
[~ 10 15 M Q ], of order supercluster size. This scale is the first to go non-linear. For 
sterile neutrinos, there is a large disparity between the two characteristic scales, so 
that perturbations with 10 13 M Q < M < 10 15 M Q are processed similarly; given an 
initial Harrison-Zel'dovich spectrum, they should all the the same final amplitude 
in linear theory. Power on scales smaller than this should be completely damped. 

In conclusion, we have proposed a candidate for warm dark matter that exists 
in the simplest extension of the standard model. Warm dark matter has several 
advantanges over cold or hot dark matter, resulting from the fact that the power 
on scales of order 1 — 5 Mpc is less than in CDM but greater than in HDM. In 
particular, the pairwise velocity dispersions in a WDM universe are likely to be 
smaller than in CDM and hence more in accord with observations. Since there is 
more power on small scales than in HDM, the epoch of galaxy formation is likely 
to be earlier and hence the observed high redshift quasars pose less of a problem 
for this model than for HDM. On large scales there is little difference among the 
three models: at present they all seem to be incompatible with the APM survey. 
Another advantage WDM has over HDM is that since the neutrino mass is higher, 
it is possible to fit more neutrinos into a given galaxy, thus evading Tremaine- 
Gunn limits 20 . Finally we point out a unique signature of WDM is an increase 
in the predicted primordial helium abundance; since a neutrino species that is in 
thermal equilibrium at the time of big bang nucleosynthesis adds AY = 0.012 to 
the primordial helium mass fraction, sterile neutrinos add 



scales smaller than the free streaming scale Xps = a Jq dt' ((p/ E) 2 } 1 ^ 2 /a(t') 



are 




(10) 



a potentially detectable deviation from the standard prediction. 
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FIGURE CAPTIONS 

1) Figure 1. Mass scales in hot dark matter and warm dark matter as a function 
of scale factor. M# (solid line) gives the mass within the horizon. Long 
dashed lines give the free streaming mass for a 30 eV (Mps,3o) and 300 
eV (Mfs^qq) neutrino. Short dashed lines are the Jeans mass for a 30 eV 
(M FS ,3o) and 300 eV (M^ Sj30 o) neutrino. 
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